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A  quantum  bit  is  a  closed  two-dimensional  Hilbert  space,  but  often  experimental  systems  have  three  or  more 
energy  levels.  In  a  Josephson  phase  qubit  the  energy  differences  between  successive  levels  differ  by  only  a  few 
percent,  and  hence  care  must  be  taken  to  isolate  the  two  desired  levels  from  the  remaining  Hilbert  space.  Here 
we  show  via  numerical  simulations  how  to  restrict  operations  to  the  qubit  subspace  of  a  three-level  Josephson 
junction  system  requiring  shorter  time  duration  and  suffering  less  error  compared  with  traditional  methods. 
This  is  achieved  by  employing  amplitude  modulated  pulses  as  well  as  carefully  designed  sequences  of  square 
wave  pulses.  We  also  show  that  tunneling  out  of  higher  lying  energy  levels  represents  a  significant  source  of 
decoherence  that  can  be  reduced  by  tuning  the  system  to  contain  four  or  more  energy  levels. 
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I.  INTRODUCTION 

The  remarkable  promise  of  quantum  computation1  has  led 
to  the  invention  of  a  significant  number  of  proposals  for 
building  a  practical  and  scalable  quantum  computer.  Several 
of  these  proposals2-6  envision  the  use  of  two  out  of  several 
energy  levels  in  a  system  as  a  quantum  bit  (qubit)  and  hence 
care  must  be  taken  to  isolate  these  levels  from  the  remaining 
Hilbert  space.  In  particular,  in  a  Josephson  phase  qubit,5'6 
energy  differences  between  successive  levels  only  differ  by  a 
few  percent.  The  challenge  in  performing  accurate  qubit  op¬ 
erations  lies  in  successfully  isolating  the  two  energy  levels 
from  the  rest  of  Hilbert  space.  In  other  words,  how  is  it 
possible  to  operate  as  quickly  and  with  as  little  error  as  pos¬ 
sible  on  the  qubit  subspace  while  isolating  the  remaining 
Hilbert  space  in  a  Josephson  phase  qubit?  This  is  especially 
important  when  the  coherence  times  of  the  system  are  short. 

A  Josephson  phase  qubit  can  be  described  by  three  energy 
levels  |0),  |l),  and  |2),  with  energies  E 0,  E  1;  and  E2,  as 
sketched  in  Fig.  1(b).  The  qubit  space  is  formed  by  |0)  and 
|l),  and  hence  we  wish  to  operate  only  within  this  subspace. 
Clearly,  the  higher-order  transition  can  be  avoided  when  ex¬ 
citing  the  a>l0  transition  by  using  a  long  enough  excitation 
duration.  However,  because  one  wants  to  maximize  the  num¬ 
ber  of  logic  operations  within  a  fixed  coherence  time,  there  is 
a  need  to  excite  the  |0)<->|l)  transition  as  quickly  as  pos¬ 
sible  without  populating  other  states. 

Here  we  numerically  model  two  techniques  which  show 
how  single-qubit  operations  can  be  improved  beyond  simple 
square  wave  pulses  (or  “hard”  pulses).  The  first  uses  ampli¬ 
tude  modulation  of  the  pulse  and  the  second  uses  composite 
pulses  that  consist  of  a  sequence  of  specially  designed  hard 
pulses.  Though  other  proposals  are  known7-9  for  implement¬ 
ing  single-qubit  rotations  within  a  subspace  of  Hilbert  space, 
this  work  takes  a  further  step  by  analyzing  decoherence  ef¬ 
fects;  specifically,  we  evaluate  the  feasibility  of  our  methods 
for  typical  parameters  of  a  Josephson  phase  qubit,  including 
the  effects  of  tunneling  out  of  higher  lying  energy  levels.  We 
show  that  these  tunneling  effects  can  be  a  significant  source 
of  decoherence  if  not  taken  into  account  properly. 


The  outline  of  the  paper  is  as  follows.  We  first  review  the 
physics  of  the  Josephson  phase  qubit  in  Sec.  II.  In  Sec.  Ill, 
we  provide  a  concrete  procedure  detailing  our  methods  to 
simulate  the  effect  of  applying  amplitude  modulated  pulses 
to  the  qubit.  Our  results  indicate  that  Gaussian  shaped 
modulation10  provides  the  best  selectivity.  In  Sec.  IV,  we 
describe  two  composite  pulses,  which  perform  better  than 
hard  pulses,  but  not  quite  as  well  as  Gaussian  shaped  pulses. 
In  Sec.  V,  we  show  how  to  include  tunneling  effects  out  of 
the  higher  energy  levels  to  estimate  the  feasibility  of  these 
techniques  in  a  real  Josephson  junction  qubit  system.  Our 
results  indicate  that  tunneling  plays  a  significant  role,  leading 
to  the  conclusion  that,  to  reduce  tunneling  effects  one  should 
use  at  least  four  energy  levels  instead  of  the  usual  three.6 

Though  our  methods  have  been  developed  in  the  context 
of  Josephson  phase  qubits,  we  believe  they  could  also  be 
fruitful  in  other  systems  where  one  wishes  to  control  a  par¬ 
ticular  subspace  of  Hilbert  space.  This  work  applies  ideas 
from  nuclear  magnetic  resonance  (NMR)  quantum 
computing,11  specifically  shaped  and  composite  pulses,  to 
other  proposals  for  implementing  a  quantum  computer;  we 
believe  this  to  be  a  very  rewarding  approach.  Furthermore, 
our  analysis  of  the  effects  of  tunneling  uses  ideas  from  quan¬ 
tum  computing  to  give  us  an  unusual  and  interesting  way  to 
understand  and  model  the  physics  of  a  real  Josephson  junc¬ 
tion  qubit. 

II.  JOSEPHSON  PHASE  QUBITS 

The  details  of  a  Josephson  phase  qubit  are  described 
elsewhere,12  and  we  review  only  the  basics  here.  The  Hamil- 
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FIG.  1.  Sketch  of  a  three  level  system  with  two  transitions  at 
frequencies  &)10  and  &)2t  ■ 
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tonian  of  the  current-biased  Josephson  junction  with  bias 
source  7,  critical  junction  current  70,  and  junction  capaci¬ 
tance  C  is 


H=2 CQ2- 


Io®0  *  I&0 

-cos  o—  — —  o. 


2tt 


2  7 T 
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where  <t>0  =  h/2e  is  the  superconducting  flux  quantum.  The 
operators  Q  and  S  correspond  to  the  charge  and  the  super¬ 
conducting  phase  difference  across  the  junction,  respectively, 
and  have  a  commutation  relationship  [S,Q]  =  2ei.  A  quan¬ 
tum  mechanical  behavior  can  be  observed  for  large  area 
junctions  in  which  70<I>0/2'7r  =  7sy8>7sc  =  e2/2C  and  when 
the  bias  current  is  slightly  smaller  than  the  critical  current 
/:£/0 .  In  this  regime  the  last  two  terms  in  77  can  be  accu¬ 
rately  approximated  by  a  cubic  potential  U(  S)  parametrized 
by  a  barrier  height  AU(I)  =  (2V2l0<f>0/3Tr)[l  —  ///0] 3/2  and 
a  quadratic  curvature  at  the  bottom  of  the  well  that  gives 
a  classical  oscillation  frequency  oop(I)  =  2l,4(2irIQ/ 
<F0C)1/2[l-7/70]1/4 

The  commutation  relation  leads  to  quantized  energy  levels 
in  the  cubic  potential.  The  quantized  energy  levels  in  this 
potential  can  be  visualized  as  indicated  by  Fig.  1(a).  Micro- 
wave  bias  currents  induce  transitions  between  levels  at  a  fre¬ 
quency  (omn  =  Emn/h  =  (Em—En)/h,  where  E„  is  the  energy 
of  state  | n).  The  two  lowest  transitions  have  frequencies 


"10  I”  T7 


5  he 


36  A  U 


"21  1-52 


10  h  co. 


36  A  U 


(2) 

(3) 


These  two  frequencies  must  be  different  to  access  the  two- 
state  system  as  a  controllable  qubit.  The  ratio  A  U/hcop  pa¬ 
rameterizes  the  anharmonicity  of  the  cubic  potential  with 
regard  to  the  qubit  states,  and  gives  an  estimate  of  the  num¬ 
ber  of  states  in  the  well. 

The  state  of  the  qubit  can  be  controlled  with  a  dc  bias 
current  Idc  and  a  time-varying  bias  current  I /JW(t)  at  fre¬ 
quency  co=(olQ,  given  by 


I(t)  =  Idc  +  Al(t)  (4) 

=  1 rfc-^w(f)cos(wt  +(/)).  (5) 

In  general,  the  Hamiltonian  for  the  lowest  three  energy 
eigenstates  of  a  Josephson  junction  system  biased  with  a 
current  7(f)  is 


The  matrix  elements  (m\S\n)  are  calculated  as  follows. 
When  1=1  dc  and  for  Idc— To,  the  system  Hamiltonian  has  a 
potential  U(S)  that  is  cubic.  We  calculate  this  Hamiltonian, 
and  solve  for  its  eigenstates  via  diagonalization,  from  which 
the  matrix  elements  (m\S\n)  can  be  simply  computed.  Cal¬ 
culating  (m|o|n)  in  this  manner,  we  obtain 
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where  the  basis  states  are  |0),  1 1),  and  |2),  from  left  to  right, 
co  is  the  frequency  of  the  applied  time-varying  current,  and 
g(t)  =  1.0147/iR,(f)  \j h/2co0lC/2  is  related  to  the  time-varying 
current  Ipw(t).  A  shaped  pulse  is  implemented  by  letting 
g(t)  vary  in  time  according  to  the  amplitude  modulation. 
The  Hamiltonian  Hm  contains  additional  diagonal  and  off- 
diagonal  elements,  but  they  are  all  sufficiently  far  enough 
off-resonance  from  w1()  and  co2\  such  that  Hm  has  only  neg¬ 
ligible  effects. 

To  calculate  the  effect  of  amplitude  modulation  it  is  con¬ 
venient  to  move  into  a  doubly  rotating  frame,  defined  by  the 
unitary  operator 
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Let  |  cf>)  =  V |  if/)  be  a  state  in  the  rotating  frame  of  V  and  |  if/)  is 
a  state  in  the  laboratory  frame.  Then  the  equation  of  motion 
for  this  state  can  be  derived  from  Schrodinger’s  equation 


ihd\  (f>)  =  H\  4>),  (10) 

where  H  is  the  rotating  frame  Hamiltonian  given  by  H 
=  V]'HV-ihV']'d,V.  This  results  in 
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where  we  have  set  E0  =  0.  In  this  work,  we  focus  only  on 
excitation  that  is  on-resonance  with  the  |0)<->|l)  transition. 
Using  El  =  hcol0,  &>=«10,  and  defining  the  energy  differ¬ 
ence  between  the  two  transitions  as  E2  —  2E  l  =  hSu),  we  ob¬ 
tain 
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We  can  now  use  this  rotating  frame  Hamiltonian  in  a 
straightforward  manner  to  numerically  calculate  the  effect 
of  hard  (square  wave)  and  shaped  pulses  on  the  three-level 
system.  We  first  discretize  the  shape  to  many  steps.  During 
each  slice  j,  we  let  the  amplitude  of  g  be  a  constant,  g  t . 
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The  unitary  evolution  in  each  slice  is  given  by  U 
=  exp(—iHjAt/h )  where  At  is  the  slice  length.  We  then  vary 
gj  from  slice  to  slice  according  to  the  modulation,  and  mul¬ 
tiply  all  unitary  evolutions  together  to  obtain  the  overall  evo¬ 
lution  of  the  system,  mathematically  described  as 


(13) 

From  U  we  can  calculate  the  leakage  out  of  the  qubit  mani¬ 
fold.  Also  note  from  Eq.  (12)  that  the  effects  of  Sol  scale  as 
the  product  tSJ2v,  the  relevant  normalized  time  unit  in  our 
problem  is  Tpw=  t SJ2tt,  and  thus  our  results  are  plotted  in 
this  way. 


III.  SHAPED  PULSES 

Shaped  pulses  are  widely  used  in  NMR  spectroscopy  and 
NMR  quantum  computing 13-1 5  because  they  can  signifi¬ 
cantly  enhance  the  selective  excitation  of  a  qubit  compared 
with  hard  pulses.  In  contrast  to  NMR,  however,  where  each 
qubit  is  represented  by  a  spin- 1/2  particle,  here  the  two  tran¬ 
sitions  share  an  energy  level,  leading  to  system  dynamics 
which  are  much  more  complex.  Although  it  is  not  immedi¬ 
ately  obvious  in  what  way  amplitude  modulating  could  be 
beneficial  in  this  system,  it  is  not  unreasonable  to  assume 
that  some  improvements  are  possible.  This  assumption  will 
be  supported  by  our  findings,  later  in  this  section. 

In  our  simulations,  we  chose  a  flip  angle  of  180°,  which 
transforms  1 0)' — > 1 1 )  and  1 1  )■ — >|0),  because  this  is  usually 
the  hardest  selective  rotation  to  achieve.  Even  though  we 
show  results  for  several  pulse  shapes,  we  only  provide  ex¬ 
plicit  functions  of  the  time  dependence  of  g(t)  for  Gaussian 
and  Hermite  180  shapes.  The  rf  envelope  g(t)  of  the  Gauss¬ 
ian  shape  is  given  by 

§ gauss(0  =  (a/tg)zx p(  -  t2/2t2g) .  (14) 

for  \t\<atg  and  g{t)  =  0  otherwise,  where  a  is  the  cut-off  of 
the  pulse  in  time  (usually  3-5),  a  is  the  amplitude  ( a 
=“1.25  for  180°  pulses  and  for  typical  values  of  a),  and  tg  is 
the  characteristic  pulse  width.  The  total  pulse  width  is  fpw 
=  2  at 

The  Hermite  180  shape  is  simply  a  Gaussian  multiplied 
by  a  second  order  polynomial.  The  RF  envelope  g(t)  of  the 
Hermite  180  shape  is 

ghrmisoW  =  (1  -  f3(t/atg)2)(a/tg)exp(  -  t2/2t2g),  (15) 

with  the  definition  of  parameters  as  before.  The  parameter  (3 
determines  how  strongly  the  Gaussian  pulse  is  modulated. 
Here,  a  =“2.2  for  a  =  3  and  1.67  for  a  =  4,  both  for  /3 
=  4  and  a  180°  flip  angle. 

Using  these  parameters  and  shapes  we  have  calculated  the 
performance  of  the  pulses.  We  define  the  error  as  e=  1 
—  1 1/ (3,3)  | 2,  where  t/ (3,3)  denotes  the  (3,3)  element  of  the 
resulting  unitary  transform.  Ideally,  £=0  if  we  are  only  op¬ 
erating  on  the  desired  qubit  subspace.  Otherwise,  the  undes¬ 
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FIG.  2.  Numerical  calculation  of  the  error  e  as  a  function  of 
normalized  pulse  width  Tvvl=tpvSJ2Tr=2atgSJ2TT  via  Fourier 
analysis  using  the  untruncated  and  truncated  Gaussian  shape,  com¬ 
pared  with  the  exact  calculation. 


ired  energy  level  is  involved  in  the  operation.  The  error  e  is 
the  same  as  the  maximum  probability  of  being  in  the  state  \2) 
after  the  application  of  the  pulse,  when  starting  from  an  ar¬ 
bitrary  superposition  ofa|0)  +  /?|l).  The  measure  e  serves  as 
a  lower  bound  error  because  in  addition  to  leakage  outside 
the  qubit  manifold,  the  desired  qubit  rotation  may  also 
slightly  deviate  from  the  ideal  rotation  even  when  e  =  0.  For 
example,  one  deviation  is  due  to  transient  Bloch-Siegert 
shifts,  similar  to  NMR,17'18  but  these  effects  can  be  corrected 
using  a  method  similar  to  the  one  described  in  Ref.  19. 

We  can  numerically  calculate  the  error  via  Eq.  (13),  but  in 
order  to  gain  an  intuitive  understanding  of  which  time- 
dependent  amplitude  modulation  could  be  most  useful,  it  is 
helpful  to  estimate  the  error  first,  using  a  more  conventional 
method  based  on  simple  bandwidth  considerations.  This  sys¬ 
tem’s  response  is  approximately  linear  for  small  rotation 
angles,  and  hence  Fourier  analysis  provides  useful  insight.  In 
fact,  Fourier  analysis  has  been  used  extensively  in  NMR 
even  for  large  rotation  angles  in  order  to  get  a  first  idea  of  the 
selectivity  of  a  shaped  pulse. 

The  relative  power  of  the  frequency  component  of  the 
untruncated  Gaussian  shaped  pulse  at  a  frequency  S0I  away 
from  oji0  is  given  by 


e(  tPw)  ~  exp[  -  ( Sj  g)2]  =  exp 


(16) 


This  is  plotted  in  Fig.  2  and  is  compared  with  both  the  Fou¬ 
rier  analysis  of  a  truncated  Gaussian  shaped  pulse  and  the 
exact  calculation  using  Eq.  (13).  For  small  pulse  widths,  the 
exact  calculation  and  the  one  based  on  the  Fourier  analysis 
of  untruncated  shapes  are  similar,  but  the  exact  calculation 
flattens  out  for  Tpw>4.  Though  approximate,  it  is  evident 
that  Fourier  analysis  still  provides  a  rough  estimate  of  the 
error,  especially  if  truncation  effects  are  included.  From  Fou¬ 
rier  analysis  we  expect  hard  pulses  to  perform  poorly  com- 
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FIG.  3.  Plot  of  the  error  e  as  a  function  of  the  normalized  pulse 
width  rpw  for  three  different  pulse  shapes  and  several  different  lev¬ 
els  of  truncation. 


pared  with  Gaussian  or  Hermite  shapes.  In  order  to  quantify 
the  performance  accurately  however,  we  must  calculate  the 
error  exactly  using  Eq.  (13). 

Figure  3  plots  e  as  a  function  of  the  normalized  pulse 
width  for  hard,  Gaussian  and  Hermite  shaped  pulses.  Clearly, 
for  low  error  rates,  a  long  pulse  must  be  used.  Since  our  goal 
is  to  apply  the  desired  rotation  as  quickly  as  possible,  it 
appears  that  the  Gaussian  shape  is  best  suited  for  this  prob¬ 
lem.  Also,  the  Gaussian  shape  could  probably  be  further  op¬ 
timized  using  our  calculational  methods  outlined  here. 

There  exist  other  pulse  shapes  which  have  been  designed 
for  NMR  experiments  to  invert  or  select  spins  over  a  very 
sharp  and  specified  bandwidth.  These  belong  to  the  class  of 
BURP  (band-selective,  uniform-response,  pure-phase) 
pulses,20  and  their  performance  is  shown  in  Fig.  4  compared 


FIG.  4.  Plot  of  the  error  e  as  a  function  of  normalized  pulse 
width  Tpw  for  several  traditional  NMR  pulse  shapes.  The  Gaussian 
shape  here  corresponds  to  <*  =  3  from  Fig.  3. 


FIG.  5.  Plot  of  the  error  e  (maximum  occupation  probability)  as 
a  function  of  time  normalized  with  respect  to  the  total  pulse  length 
tpw  during  a  hard  and  Gaussian  shaped  pulse  and  the  second  com¬ 
posite  pulse  (Sec.  IV  B).  A  Gaussian-shaped  pulse  with  a  pulse 
width  of  rpw=4  has  a  very  small  error  (about  1CT8),  but  during  the 
pulse  the  error  can  be  as  high  as  5X  lCT2.  The  error  during  the 
pulse  can  be  reduced  by  applying  the  pulse  for  a  longer  duration 
since  the  maximum  error  scales  as  approximately  l/rpw ,  based  on  a 
Fourier  analysis  of  the  truncated  pulses. 

with  hard,  Gaussian  and  Hermite  shapes.  Rather  surprisingly, 
these  specially  designed  shapes  perform  rather  poorly  in  our 
three-level  system.  This  is  probably  because  these  pulses 
were  designed  for  use  in  two-level  systems  whose  dynamics 
are  quite  different  compared  to  a  Josephson  junction  where 
successive  transitions  share  one  energy  level. 

It  is  also  interesting  to  note  that  the  error  e  (or  the  maxi¬ 
mum  occupation  probability  of  state  |2))  during  the  applica¬ 
tion  of  the  pulse  may  be  much  higher  than  at  the  end  of  the 
pulse,  as  indicated  in  Fig.  5.  Depending  on  the  experiment, 
this  may  still  be  undesirable,  for  example  when  the  lifetime 
of  the  third  energy  level  is  short.  We  discuss  the  impact  of 
such  short  lifetimes  later,  in  Sec.  V. 

While  pulse  shaping  clearly  offers  advantages  in  this 
three-level  system,  it  may  not  necessarily  be  straightforward 
to  do  so  experimentally,  and  hence  it  is  useful  to  look  at 
alternatives  to  pulse  shaping.  In  the  next  section  we  discuss 
how  composite  pulses  could  provide  an  alternative  approach 
towards  improving  single -qubit  operations  in  a  Josephson 
phase  qubit. 

IV.  COMPOSITE  PULSES 

A  composite  pulse  is  a  pulse  that  consists  of  a  sequence  of 
individual  pulses,  and  is  designed  to  reduce  certain  types  of 
errors,  but  often  at  the  cost  of  longer  total  duration.  The 
individual  pulses  are  typically  square  wave  pulses,  but  can  in 
principle  be  shaped  pulses  as  well.  Similar  to  shaped  pulses, 
composite  pulses  have  also  found  a  wide  variety  of  applica¬ 
tions  in  NMR.13’21'22  The  overall  unitary  evolution  U  of  a 
composite  pulse  is  calculated  via  U=YVkUk,  where  the  Uk 
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are  the  unitary  evolutions  of  the  individual  pulses  contained 
in  the  sequence.  We  next  describe  the  design  of  two  different 
composite  pulses,  both  of  which  perform  better  than  hard 
pulses. 


A.  Composite  pulses — method  1 

We  begin  by  noting  that  g  =  0  is  a  valid  choice  for  a  pulse. 
This  corresponds  to  applying  no  electromagnetic  radiation 
for  some  time.  Working  with  the  rotating  frame  Hamiltonian 
of  Eq.  (12),  and  letting  the  quantum  system  freely  evolve  for 
a  time  t  =  tt!  8w  ,  we  obtain  the  transformation 


Z2  = 


1  0 
0  1 
0  0 


0 

0 

-1 


(17) 


where  the  subscript  denotes  that  state  \2)  acquires  a  180° 
phase  shift  with  respect  to  all  other  states.  Now,  if  we  sand¬ 
wich  a  driven  evolution  U{t)  between  two  free  evolution 
periods  Z2,  we  obtain 


/  . 

It 

~h 

'0 

g 

0 

\ 

V(t)  =  Z2U(t)Z2  =  ex  p 

8 

0 

-V2g 

\ 

.0 

-V2g 

h  8U>  . 

/ 

(18) 

where  we  assumed  g  to  be  real  and  U(t)=e  ,  Note 
that  we  can  interpret  the  exponent  of  V(t)  as  applying  an 
x  rotation  on  the  |0)<->|l)  transition  and  a  negative 
x-rotation  on  the  |l)<->|2)  transition,  in  the  limit  of 
fast  pulses  (8^/g  — >0).  If  we  now  apply  a  second  pulse 
U(t)  which  can  be  regarded  as  an  x-rotation  on  both  transi¬ 
tions,  the  excitation  of  the  undesired  transition  is  undone. 
The  overall  transformation  is  then  given  by  W(t ) 
=  U(t/2)V(t/2)  =  U(t/2)Z2U(t/2)Z2,  and  is  an  interesting 
pulse  sequence  because  it  becomes  a  nontrivial  unitary  trans¬ 
form  on  just  the  desired  qubit  manifold  in  the  limit  of  fast 
pulses  with 

This  procedure,  however,  only  works  for  small  rotation 
angles.  Even  though  the  interpretation  of  the  exponents  of 
V(t)  and  U(t)  above  is  correct,  the  matrices  V{t),  U{t),  and 
W[t )  themselves  have  elements  connecting  the  |0)  and  \2) 
states  whose  magnitudes  are  second  order  in  time  t.  We  can 
only  ignore  these  for  small  rotation  angles.  For  larger  rota¬ 
tion  angles,  one  may  have  to  apply  a  sequence  of  W{t/n), 


R(t)  =  [W{t/n)]n,  (19) 

to  suppress  the  |0)<->|2)  excitation.  Furthermore,  if  we 
chose  n  sufficiently  large,  then  state  |2)  also  only  has  a  very 
small  transient  population  during  the  pulse  sequence.  The 
total  operation  time  is  equal  to  n(2  tt/  S^)  +  t  (where  ris  the 
total  duration  of  the  electromagnetic  radiation),  which  be¬ 
comes  quite  long  for  large  n.  Surprisingly,  this  composite 
method  still  performs  rather  well  even  when  the  system  is  far 
away  from  the  short  pulse  limit  g  S>  Su ,  as  shown  in  Fig.  6. 
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FIG.  6.  Plot  of  the  error  e  as  a  function  of  normalized  total  time 
rpw  for  the  two  composite  pulses  (n  =  2)  described  in  Sec.  IV,  com¬ 
pared  with  a  hard  and  Gaussian  shaped  pulse.  The  first  composite 
method  outperforms  hard  pulses  whenever  the  pulses  are  multiples 
of  2tt/ Su  even  though  the  limit  of  does  not  apply  here.  For 
the  second  composite  pulse,  when  the  applied  pulses  are  integer 
multiples  of  2  tt/  <5w  ,  the  error  is  minimized  as  expected. 


B.  Composite  pulses — method  2 

Another  composite  pulse  design  is  based  on  the  previous 
solution,  supplemented  by  Bloch-sphere  intuition  and  a 
knowledge  of  the  excitation  profile  of  a  hard  pulse.  Suppose 
we  excite  a  two-level  system  SJ2ir  hertz  away  from  reso¬ 
nance  via  a  hard  pulse  of  duration  r.  Whenever  8J2tt  is  an 
integer  multiple  of  1/r,  the  system  is  not  excited.  However, 
in  our  system  the  two  transitions  share  an  energy  level. 
Hence,  even  if  we  excite  the  undesired  transition  off- 
resonance,  there  is  still  a  substantial  error,  as  evident  from 
Fig.  3.  Nonetheless,  such  a  carefully  timed  pulse  is  still  use¬ 
ful  for  the  design  of  our  composite  pulse. 

For  the  first  step  of  the  composite  pulse,  let  the  hard  pulse 
be  applied  on  resonance  with  the  1 0)+— >  1 1 )  transition  for  a 
time  of  2tt/  Su.  Let  this  pulse  be  denoted  by  U (2  it/  S0>)  and 
let  the  power  be  such  that  we  would  obtain  only  a  small 
rotation  angle  (less  than  45°).  The  resulting  matrix  element 
of  U{2tt/ <5W)  connecting  the  |l)  and  \2)  states  is  very  small 
because  this  transition  is  SJ2tt  Hz  off-resonance,  in  accor¬ 
dance  with  our  Fourier  intuition  from  the  previous  para¬ 
graph.  However,  U (2  tt/  80l)  still  has  a  matrix  element  con¬ 
necting  the  |0)  and  |2)  states  because  the  two  transitions  share 
an  energy  level,  leading  to  errors  (similar  to  U{t)  and  V(  t) 
from  above).  The  last  step  of  the  composite  pulse  simply 
consists  of  undoing  the  |0)<->|2)  excitation. 

The  |0)  <->|2)  excitation  can  be  reversed  by  using  the  fact 
that  a  180°  z  rotation  sandwiched  between  two  x -rotations 
leads  to  no  net  excitation.  Suppose  we  apply  the  sequence 
C=  U(2  tt/ S0>)Z2U(2  tt/S^).  Let  us  investigate  the  effects 
of  C  on  the  three  possible  excitations:  (1)  The  first  and  last 
pulses,  U{2tt/ <5W),  excite  the  |0)<->|l)  transition,  whereas 
Z2  has  no  effect  on  the  states  |0)  and  |l).  (2)  Neither 
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U(2tt/  S a)  nor  Z2  excite  the  1 1 ) <->  1 2 )  transition  as  described 
in  the  previous  paragraph.  (3)  The  two  pulses  U{2tt/S1i1) 
lead  to  excitations  between  the  states  |0)  and  |2),  which  can 
be  regarded  as  x -rotations,  whereas  Z2  acts  like  a  180° 
z-rotation  on  the  same  states.  Hence,  the  sequence  C  leads  to 
no  net  excitation  between  the  |0)  and  \2)  states. 

As  a  result,  the  sequence  C  acts  non-trivially  only  on  the 
desired  qubit  manifold.  By  carefully  adjusting  the  power  of 
the  pulses  we  can  create  any  arbitrary  x  rotation  on  the  de¬ 
sired  transition.  The  sequence  thus  consists  of  the  application 
of  two  pulses  of  length  2  tt!  Sio  interlaced  with  a  free  evolu¬ 
tion  of  length  7 t/  Sol . 

Similar  to  the  previous  composite  pulse,  this  method  re¬ 
sults  in  low  errors  e  only  when  U (2  tt/  S(j  implements  small 
rotations.  Hence  in  general,  we  may  also  have  to  apply  n 
instances  of  U(2tt/  Su1)Z2U(2tt/  Sw)  if  our  goal  is  to  achieve 
large  rotation  angles.  The  resulting  duration  of  this  compos¬ 
ite  pulse  is  2.5«(27t/<5w). 

Figure  6  shows  the  results  of  the  composite  pulse  methods 
compared  to  a  hard  pulse  and  the  Gaussian  shaped  pulse.  The 
time  axis  denotes  the  total  duration  of  the  pulses.  For  the 
composite  pulses,  this  includes  the  duration  of  the  electro¬ 
magnetic  radiation  plus  the  delay  period  implementing  Z2. 
As  can  be  seen,  the  two  composite  pulse  methods  outperform 
hard  pulses  but  not  a  Gaussian  pulse  with  a=  3.  It  is  clear 
that  the  two  simple  composite  pulses  may  provide  useful 
alternatives  for  implementing  accurate  single  qubit  rotations 
in  a  Josephson  phase  qubit.  Furthermore,  it  should  be  pos¬ 
sible  to  combine  shaped  and  composite  pulses  to  even  further 
improve  qubit  operations. 

Nonetheless,  as  is  the  case  with  shaped  pulses,  the  occu¬ 
pation  probability  of  state  \2)  during  the  application  of  the 
composite  pulse  can  be  rather  large.  For  the  second  compos¬ 
ite  pulse,  this  number  is  3  X  10~2  (see  Fig.  5),  which  is  much 
larger  than  the  occupation  probability  at  the  end  of  the  pulse. 
In  the  next  section  we  show  how  transient  populations  in  the 
third  energy  level  can  be  highly  undesirable  in  the  presence 
of  high  tunneling  rates. 

V.  EFFECTS  OF  TUNNELING 

Thus  far  we  have  only  considered  the  ideal  case  where  the 
two  transitions  of  the  three-level  system  are  close  in  fre¬ 
quency  but  otherwise  there  exist  no  sources  of  decoherence. 
However,  in  a  real  Josephson  phase  qubit,  the  quantum  state 
can  tunnel  through  the  barrier  of  the  cubic  well,  and  this 
process  acts  as  a  source  of  decoherence.  What  consequences 
do  we  expect  from  this? 

From  first  principles  we  know  that  the  tunneling  rate  de¬ 
pends  exponentially  on  the  barrier  height  and  width.  Hence 
in  our  cubic  potential,  we  expect  the  upper  level  to  be  most 
susceptible  to  tunneling.  When  the  tunneling  rate  out  of  this 
energy  level  is  high,  then  there  exists  a  significant  source  of 
decoherence  if  the  state  becomes  even  transiently  populated. 
It  is  thus  desirable  to  keep  transient  populations  in  state  \2)  as 
small  as  possible  during  our  single  qubit  rotations.  This  may 
not  necessarily  be  the  case  however,  as  indicated  in  Fig.  5. 
We  next  estimate  the  impact  of  this  tunneling  effect,  first  by 
using  a  simple  model,  followed  by  a  more  rigorous  approach. 


FIG.  7.  Plot  of  P,  as  a  function  of  normalized  pulse  width  for  a 
hard,  Gaussian  and  composite  pulse  taking  T SJ2tt 
^(a>10/2ir)/100.  Clearly,  the  estimate  of  the  tunneling  probability  is 
several  orders  of  magnitude  higher  than  e,  indicating  that  tunneling 
effects  are  important. 


A.  A  simple  tunneling  model 

A  first  model  that  gives  insight  into  the  importance  of 
transient  populations  during  the  pulses  can  be  set  up  as  fol¬ 
lows.  Let  the  probability  of  being  in  state  \2),  and  the  tun¬ 
neling  rate  out  of  \2)  be  defined  as  p2  and  F 2  respectively. 
The  probability  of  tunneling  out  of  state  |2)  can  then  be 
calculated  by  P,  =  Jp2T2dt.  If  we  bias  the  system  such  that 
r2  is  on  the  order  of  SJ2tt,  which  is  about  10  times  larger 
than  the  inverse  of  typical  pulse  widths,6  we  expect  any  tran¬ 
sient  populations  in  state  \2)  to  tunnel  out  of  the  potential 
well  during  the  pulse.  Since  the  tunneling  rates  out  of  1 1 ) 
(T !)  and  |0)  (T 0)  are  about  103  and  106  times  less  than  T2, 
we  ignore  their  effects  in  this  simple  calculation.6 

In  Fig.  7  we  plot  P,  as  a  function  of  normalized  pulse 
width  for  a  hard,  Gaussian,  and  the  second  composite  pulse 
with  a  tunneling  rate  r2'=5(&)10/27r)/100®s SJ2tt.  Clearly, 
the  error  is  several  orders  of  magnitude  higher  than  e,  and 
tunneling  thus  appears  to  be  an  important  source  of  decoher¬ 
ence  for  this  system,  as  expected.  Note  that  P,  is  only  an 
approximate  overall  error  because  it  reflects  only  the  prob¬ 
ability  of  tunneling  and  does  not  include  the  occupation 
probability  of  being  in  state  |2)  at  the  end  of  the  pulse,  which 
can  be  non-zero  regardless  of  the  tunneling  rate.  We  next 
describe  a  more  rigorous  approach  that  analyzes  the  effects 
of  tunneling  and  includes  the  leftover  occupation  probability 
of  state  |2). 

B.  A  tunneling  model  using  the  operator-sum  representation 

We  shall  model  the  tunneling  behavior  similar  to  ampli¬ 
tude  damping  by  using  the  operator-sum  representation.23 
This  type  of  model  has  been  successfully  used  to  predict  the 
impact  of  decoherence  in  several  NMR  quantum  computing 
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experiments.15'24  In  the  operator  sum  representation,  an  ini¬ 
tial  density  matrix  p,  is  mapped  to  a  final  density  matrix  (u 
via 


P/=2  Ekp,El  (20) 

k 


where  H,kE\Ek=I  and  Ek  are  the  Kraus  operators.23  For  am¬ 
plitude  damping  acting  on  a  single  qubit,  we  have  only  two 
Kraus  operators,  taking  the  form 


1  0 

°~  0  sjl-PT 

'0  #7 


(21) 


(22) 


where  P r  =  1  —  e~,r  with  I  denoting  the  inverse  lifetime  of 
the  |l)  state.  From  this,  it  becomes  clear  that  any  quantum 
state  eventually  collapses  to  the  ground  state  |0). 

The  tunneling  mechanism  out  of  state  |2)  is  now  modeled 
by  using  a  fourth,  fictional  auxiliary  energy  level  |7)  which 
acts  as  a  reservoir  for  the  tunnel  states.  We  now  define  a 
fictitious  qubit  with  basis  states  |2)  and  |7).  The  tunneling 
mechanism  from  |2)  to  the  auxiliary  level  is  then  captured  by 
modeling  amplitude  damping  on  the  fictitious  qubit.  In  this 
case  the  two  Kraus  operators  take  the  forms 

‘  1 

0 


0 


0  0  0 

1  0  0 

0  \]l-Pr  0 

0  0  1 


(23) 


E  i  = 


0  0 
0  0 
0  0 
0  0 


0  O' 
0  0 
0  0 
Vfv  o_ 


(24) 


with  the  basis  states  |0),  |l),  |2),  and  |T)  going  from  left  to 
right.  These  operators  also  ensure  that  if  the  state  tunneled, 
then  all  coherences  to  that  state  vanish.  For  example,  if  the 
initial  state  was  an  equal  superposition  of  |0)  and  |2),  then 
after  tunneling  we  have  a  mixed  state  without  coherences, 
and  find  that  the  system  is  in  state  |0)  or  |T),  each  with 
probability  0.5. 

The  tunneling  rates  F  ,  and  I  ()  are  ~  103  and  ~  106  times 
lower  than  T2,  as  described  earlier.  When  Y ^  SJItt,  we 
estimate  their  contributions  to  be  only  on  the  order  of  1 0  5 
for  typical  pulse  widths,  using  the  simple  tunneling  model 
from  above.  Hence  we  do  not  include  their  effects  in  our 
operator-sum  approach.  Later,  we  show  how  to  generalize 
our  method  to  include  tunneling  from  an  arbitrary  number  of 
levels,  and  use  this  generalized  method  to  verify  that  Y ,  and 
T0  are  indeed  negligible,  though  we  do  not  explicitly  show 
our  results  here. 

The  Kraus  operators  above  describe  the  tunneling  mecha¬ 
nism,  but  they  do  not  include  excitation  effects.  The  next 


FIG.  8.  Plot  of  the  error  as  a  function  of  normalized  pulse  width 
when  including  tunneling  effects  with  r2  =  <5w*s(&)1o/2Tr)/100.  The 
error  is  given  by  the  maximum  of  |2)(2|  +  |T)(T|  over  all  possible 
input  superpositions  a|0)  +  fe|  1 ).  Now,  the  hard  pulse  shows  the 
best  performance,  but  still  with  a  large  error. 


step  to  model  the  single-qubit  rotations  in  the  presence  of 
tunneling  is  to  include  the  excitation  of  the  transitions  due  to 
the  applied  radiation.  The  excitations  of  the  two  transitions  in 
our  three-level  system  do  not  commute  with  the  tunneling 
mechanism.  Thus,  strictly  speaking,  one  may  have  to  derive 
a  more  complex  formalism  that  captures  both  the  tunneling 
effect  and  the  qubit  excitation  simultaneously.  However,  we 
approximate  the  simultaneous  behavior  by  slicing  the  micro- 
wave  pulses  into  many  discrete  steps  (as  has  been  done  in 
Secs.  Ill  and  IV),  and  then  simulate  the  tunneling  and  exci¬ 
tation  one  after  the  other  in  each  step.  In  the  limit  as  the 
number  of  steps  goes  to  infinity  this  approximation  becomes 
exact;  we  have  used  256  steps  for  our  simulation. 

From  Eq.  (20),  it  becomes  clear  that  the  operator  formal¬ 
ism  takes  density  matrices  as  the  input,  which  we  have  not 
yet  specified.  Our  goal  is  to  find  the  maximum  possible  error 
over  all  possible  input  states,  which  are  arbitrary  superposi¬ 
tions  of  the  qubit,  defined  as  a\0)  +  b\  1 ).  We  vary  the  values 
a  and  b  to  maximize  the  error,  which  is  defined  as  the  sum  of 
the  |2)(2|  and  |  T)(T |  elements  of  the  resulting  density  matrix. 
This  corresponds  to  the  maximum  probability  of  being  in  the 
states  \2)  and  \T).  It  is  a  useful  error  measure,  and  in  the 
absence  of  tunneling  it  yields  the  same  error  e  that  we  have 
defined  before. 

In  Fig.  8  we  plot  the  error  as  a  function  of  normalized 
pulse  width  for  a  tunneling  rate  F2SS5  Wk/IOO®3  SJItt.  From 
this  figure,  we  notice  that  the  error  is  on  the  order  of  1CT2  for 
typical  pulse  widths.  Note  that  especially  for  the  Gaussian 
shaped  pulse,  the  exact  calculation  matches  the  result  of  the 
simple  tunneling  model  from  Fig.  7  reasonably  well.  It  is 
evident  from  both  plots  that  for  typical  pulse  widths  the  error 
is  much  larger  than  the  10  4  threshold  required  for  fault- 
tolerant  computations.25  How  can  we  reduce  the  error  result¬ 
ing  from  tunneling  without  significantly  increasing  the  pulse 
widths? 
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From  the  simple  theoretical  model  described  earlier,  the 
most  straight  forward  method  is  to  decrease  the  tunneling 
rates,  which  can  be  done  by  adjusting  the  bias  current  I  to 
include  more  than  three  energy  levels  in  the  well.  In  this 
case,  the  frequency  difference  8J2jt  between  the  transitions 
becomes  slightly  smaller.  However,  the  tunneling  rate  is  ex¬ 
ponentially  dependent  on  the  barrier  thickness,  and  thus  we 
can  reduce  the  tunneling  rate  by  many  orders  of  magnitude 
while  only  slightly  decreasing  the  frequency  difference  be¬ 
tween  the  transitions.  With  four  levels  in  the  well,  we  can 
reduce  V  2  by  three  orders  of  magnitude  to  about 
(  SJ2it)/1000  while  decreasing  8J2tt  by  only  30%.  Though 
r3  is  still  high  (about  8J2jt),  level  |3)  is  sparsely  populated 
because  the  |2)<->|3)  transition  is  about  28J2rnWx  off- 
resonance  and  hence  should  not  be  significantly  excited. 

The  tunneling  effects  from  \2)  and  |3)  during  single-qubit 
rotations  can  be  calculated  via  the  operator-sum  representa¬ 
tion  as  follows.  The  effect  of  the  radiation  and  tunneling 
acting  simultaneously  is  approximated  as  before  by  slicing 
the  pulse  into  many  steps  and  simulating  tunneling  and  ex¬ 
citation  one  after  the  other  in  each  step.  The  tunneling  is 
modeled  by  sequentially  modeling  tunneling  out  of  states  |2) 
and  |3).  In  other  words,  we  first  apply  the  two  Kraus  opera¬ 
tors  that  describe  the  tunneling  from  \2)  to  \T).  Then,  we 
apply  the  two  Kraus  operators  that  the  describe  the  tunneling 
from  1 3)  to  |7).  Strictly  speaking  this  is  not  correct  because 
the  Kraus  operators  that  describe  tunneling  from  \2)  do  not 
commute  with  those  describing  the  tunneling  from  |3).  How¬ 
ever,  since  the  pulse  is  already  discretized  into  many  steps 
(we  used  256),  during  each  of  which  we  model  tunneling 
separately,  our  approach  becomes  a  good  approximation.  The 
resulting  four  Kraus  operators  can  be  easily  derived  from 
Eqs.  (23)  and  (24).  Furthermore,  we  can  model  the  tunneling 
from  an  arbitrary  number  of  states  in  this  manner,  and  have 
verified  that  the  error  for  a  three -level  system  is  indeed  domi¬ 
nated  by  the  tunneling  from  state  \2)  and  that  tunneling  from 
|l)  and  |0)  is  negligible. 

Figure  9  shows  the  error  for  a  four  level  system  with 
tunneling  rates  r2s“( <5^/2 7r)/ 1000,  and  r3«=  8J2tt.  As  can 
be  seen,  the  results  are  much  improved  compared  with  Fig. 
8.  In  fact,  now  the  composite  pulse  can  be  as  good  as  the 
Gaussian  shaped  pulse. 

It  is  possible  to  continue  increasing  the  number  of  levels 
in  the  well  to  suppress  tunneling  even  further.  However,  this 
may  not  be  practical  much  beyond  four  or  five  levels  because 
the  state  measurement  needs  a  transition  between  state  |l) 
and  a  higher  energy  state  with  a  very  large  tunneling  rate.6 
This  becomes  increasingly  difficult  for  more  energy  levels  in 
the  well.  In  the  four  level  case,  it  is  possible  to  directly  excite 
the  1 1  )•<— >  1 3 )  transition.26 

Finally,  we  point  out  that  other  proposals  for  performing 
single  qubit  rotations  in  a  three  or  multilevel  system7  9  are 
also  problematic  whenever  energy  levels  with  high  tunneling 
rates  are  populated.  We  believe  that  our  method  for  simulat¬ 
ing  the  tunneling  mechanism  could  be  useful  to  estimate  the 
feasibility  of  these  and  other  methods. 

VI.  CONCLUSIONS 

In  summary,  we  have  shown  how  shaped  and  composite 
pulses  can  improve  the  accuracy  of  single-qubit  operations  in 


FIG.  9.  Plot  of  the  error  as  a  function  of  normalized  pulse  width 
when  including  tunneling  effects  in  a  four  level  system.  The  tunnel¬ 
ing  rate  T2  is  (  8J2tt)/\  000,  and  T3  is  1000  times  higher.  The  error 
is  given  by  |2)(2|  +  |3)(3|  +  |7')(7’|  after  maximizing  over  all  pos¬ 
sible  input  superpositions  a|0)  +  fo|l).  In  this  case,  the  composite 
pulse  can  perform  about  as  well  as  a  Gaussian  pulse. 

a  Josephson  phase  qubit.  We  estimate  the  feasibility  of  our 
methods  by  including  tunneling  effects  to  show  that  tunnel¬ 
ing  can  be  a  significant  source  of  decoherence,  and  we  con¬ 
clude  that  operating  a  Josephson  phase  qubit  with  three  lev¬ 
els  is  not  recommended.  Instead,  to  reduce  decoherence 
effects  from  unwanted  tunneling  one  may  wish  to  operate  the 
system  with  four  energy  levels.  Note  that  we  have  ignored  all 
other  sources  of  decoherence  in  our  calculations  such  as  de¬ 
coherence  due  to  current  noise  and  dissipation12'27  which 
would  act  in  addition  to  leakage  and  tunneling. 

Our  proposed  method  of  using  shaped  and  composite 
pulses  should  be  realizable  using  commercial  electronics. 
Typical  junction  parameters  are  10  GHz,  A U/hcop 

=  4,  and  leading  to  typical  pulse  widths  fpw 

that  are  on  the  order  of  a  few  nanoseconds  or  tens  of  nano¬ 
seconds.  Using  analog  filters,  commercial  electronics  can  ac¬ 
commodate  shaped  and  composite  pulses  at  even  shorter 
time  scales  than  nanoseconds.  Considering  that  reported  co¬ 
herence  times  are  100  ns  or  longer,5'26  the  performance  of 
shaped  and  composite  pulses  can  be  immediately  experimen¬ 
tally  investigated. 

In  this  work  we  have  only  considered  single-qubit  gates 
which  are  not  sufficient  to  build  a  quantum  computer.  Hence, 
to  estimate  the  feasibility  of  building  a  quantum  computer 
using  coupled  Josephson  junctions,  it  will  be  necessary  to 
characterize  the  performance  of  two-qubit  gates  as  well. 
Work  towards  that  direction  has  already  begun,28  but  without 
considering  pulse  shaping.  As  we  have  demonstrated,  shaped 
pulses  are  capable  of  reducing  leakage  to  other  energy  levels 
and  hence  we  envision  that  our  techniques  can  also  offer 
practical  solutions  for  two-qubit  gates. 

Our  work  demonstrates  that  NMR  quantum  computing 
and  spectroscopy  techniques  are  a  valuable  resource  when 
applied  to  related  quantum  systems,  for  example  the  super¬ 
conductor  based  qubit  system.  We  believe  that  continued 
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effort  towards  this  direction  could  prove  fruitful  for  other 
implementations  of  quantum  computers  including  solid  state 
and  trapped  ion  implementations.  Furthermore,  we  applied 
fundamental  ideas  from  quantum  computing  to  simulate  tun¬ 
neling  effects  in  Josephson  junction  qubits,  illustrating  how 
quantum  computing  is  useful  in  modeling  and  simulating  the 
physics  of  real  quantum  systems  in  a  straightforward  manner. 
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